UM-TH-00-30 
MZ-TH/00-33 
CLNS-00/1686 



hep-ph/0008T56| 



Top quark mass definition and 
tt production near threshold at the NLC 

Oleg Yakovlev^'" and Stefan Groote^'^ '' 

^ Randall Laboratory of Physics, University of Michigan, 
Ann Arbor, Michigan 48109-1120, USA 

^ Institut fiir Physik der Johannes-Gutenberg- Universitat, 
Staudinger Weg 7, 55099 Mainz, Germany 

^ Floyd R. Newman Laboratory, Cornell University, 
Ithaca, NY 14853, USA 



Abstract 

We suggest an infrared-insensitive quark mass, defined by subtracting the soft part of the 
quark self energy from the pole mass. We demonstrate the deep relation of this definition 
with the static quark-antiquark potential. At leading order in 1/m this mass coincides 
with the PS mass which is defined in a completely different manner. Going beyond static 
limit, the small normalization point introduces recoil corrections which are calculated here 
as well. Using this mass concept and other concepts for the quark mass we calculate the 
cross section of e^e~ ti near threshold at NNLO accuracy adopting three alternative 
approaches, namely (1) fixing the pole mass, (2) fixing the PS mass, and (3) fixing the 
new mass which we call the PS mass. We demonstrate that perturbative predictions 
for the cross section become much more stable if we use the PS or the PS mass for the 
calculations. A careful analysis suggests that the top quark mass can be extracted from 
a threshold scan at NLC with an accuracy of about 100 — 200 MeV. 
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1 Introduction 



One of the main goals of future e~^e~ and fi^fi^ colliders such as the Next Linear Collider (NLC) 
and the Future Muon Collider (FMC) will be to measure and to determine the properties of the 
top quark which was first discovered at the Tevatron with a mass of m = 174.3 ± 5 GeV [|2|. 
Although the top quark will be studied at the LHC and the Tevatron (RUN-II) with an expected 
accuracy for the mass of 2— 3 GeV, the most accurate measurement of the mass with an accuracy 
of 0.1% (100 - 200 MeV) is expected to be obtained only at the NLC §. 

Due to the large top quark width, the top-antitop pair cannot hadronize into toponium 
resonances. The cross section appears therefore to have a smooth line-shape showing only a 
moderate IS" peak. In addition the top quark width serves as an infrared cutoff |^ and as 
a natural smearing over the energy As a result, the nonperturbative QCD effects induced 
by the gluon condensate are small |]^, |^, allowing us to calculate the cross section with high 
accuracy by using perturbative QCD even in the threshold region. Many theoretical studies 
at LO and NLO have been done in the past for the total cross section ^, ^ ^, for the 
momentum distribution [0, |ll], also accounting for electro- weak corrections [|l^, |13|, [l^, and 



for the complete NLO correction including non-factorizable corrections |T^, IB, 17, TE\. 

Recently, the NNLO analysis of the inclusive threshold production cross section has been 
reported [| 



^0 21 



|25[| . The results of the NNLO analysis are summarized in a 
review article [^. To summarize the results for a standard approach using the pole mass, 
the NNLO corrections are uncomfortably large, spoiling the possibility for the top quark mass 
extraction at NLC with good accuracy because the IS" peak is shifted by about 0.5 GeV by the 
NNLO, the last known correction. One of the main reasons for this is the usage of the pole 
mass in the calculations. It was realized that such type of instability is caused by the fact that 
the pole mass is a badly defined object within full QCD [^, |28|| . 

In this paper we suggest a definition of the quark mass alternative to the pole mass. We 
call it PS mass because in the static limit this mass coincides with the potential subtracted 
(PS) mass, even though the PS mass is defined in a different manner This is the reason 
why we include both of them in this paper. The ratio of the small normalization point and the 
mass introduce relativistic corrections (called recoil corrections) which result in the difference 
between the two masses at higher orders in 1/m. In contrast to the pole mass, the PS and 
the PS mass are not sensitive to non-perturbative QCD effects. We derive recoil corrections 
to the relation of the pole mass to the PS mass and demonstrate that perturbative predictions 
for the cross section become much more stable at higher orders of QCD (shifts are below 
0.1 GeV) if we use the PS or PS mass for the calculations, as it is the case for any of the other 
threshold masses. This understanding removes one of the obstacles for the accurate top quark 
mass measurement and it can be expected that the top quark mass will be extracted from a 
threshold scan at NLC with an accuracy of about 100 — 200 MeV. The necessity of isolating the 
IR contributions in the mass calculation and the consideration of a short distance mass have 
been studied intensively ^ ^ |3^. The applications of the PS, LS and 15* mass have 

and reviewed in Ref. 



been reported recently by several groups 
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26 . The 



detailed comparison of our results using the PS mass with results of other groups have been 
performed in Ref. ^6] so that we refer the reader to this reference for details. 

This paper is therefore devoted to two subjects. First, we define the PS mass which allows 
us to calculate recoil corrections near the threshold. Second, we calculate the cross section 
for e^e~ tt near threshold with NNLO accuracy using three alternative mass concepts, i.e. 
(1) the pole mass, (2) the static PS mass, and (3) the new PS mass. The results of the first 
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approach has aheady been reported in Ref. ||2T|, the preliminary resuhs for the top quark pair 



production near threshold using the static PS mass were presented at the workshop "Physics at 
Linear CoUiders" ^ by one of the authors. The results on the second and third alternative 
concepts which we present here are new. 

The paper is organized as follows: In Sec. 2 we give the definition of the PS mass in terms 
of the soft part of the self energy and calculate the leading order contribution. In Sec. 3 we 
calculate the two- loop corrections and present our final result. In Sec. 4 we use the obtained 
results for the analysis of the top quark pair production. In Sec. 5 we give our conclusions. 
The Appendix contains explicit expressions for the coefficients which occur in the text. 



2 Mass definitions 



The top quark mass is an input parameter of the Standard Model. Although it is widely 
accepted that the quark masses are generated due to the Higgs mechanism, the value of the mass 
cannot be calculated from the Standard Model. Instead, quark masses have to be determined 
from the comparison of theoretical predictions and experimental data. 

It is important to stress that there is no unique definition of the quark mass. Because 
the quark cannot be observed as a free particle like the electron, the quark mass is a purely 
theoretical notion and depends on the concept adopted for its definition. The best known 
definitions are the pole mass and the MS mass. However, both definitions are not adequate 
for the analysis of top quark production near threshold. The pole mass should not be used 
because it has the renormalon ambiguity and cannot be determined more accurately than 
400 MeV P, [ 



300 



(see also Refs. p9 



[]). The MS mass is an Euclidean mass, defined at 
high virtuality, and therefore destroys the non-relativistic expansion. Instead, it was recently 
suggested to use threshold masses like the low scale (LS) mass [|^, the potential subtracted 
(PS) mass P3, or one half of the perturbative mass of a fictions 1^5*1 ground state (called IS" 



mass) p3[. In this paper we study the static PS mass suggested in Ref. [p^] 



mps 



m 



pole 



— 6mps with 6m 



PS 



Vcm) 



where Vc is the quark- ant iquark Coulomb potential. In order to understand why this mass 
definition is more adequate than the pole mass and to see that the pole mass is very sensitive 
to long distance effects, it is enough to consider the one-loop expression for the self energy 
diagram. Taking the residue in ko, one obtains a soft self energy contribution which comes 
from momenta k with \k\ < Aqcd, 



5m = AnasCp J 



d^k 



(27r)3 2|^|2 

fc|<AQCD 



QCD- 



(2) 



We observe that the pole mass has a non-perturbative uncertainty of order Aqcd which then 
penetrates into consequent perturbative QCD calculations. It is easy to realize that the PS 
mass is free of this ambiguity. Indeed, the term 6m in Eq. (H) cancels in the definition of the 
PS mass as given in Eq. (I]) as well as in the combination 2mpoie + V{r). The definition in 
Eq. (0) has been given in Ref. |29| but has already been discussed implicitly in Ref. |28]. The 
remarkable step made in Ref. ||2^ is to use this definition beyond one-loop order. It has been 
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Figure 1: The general structure of the self energy diagram of a quark 



proven in Ref. that the cancellation of the infrared QCD contributions to the PS mass in 
Eq. (|l|) holds even at higher loop orders. 

In this paper we suggest a natural definition of a short distance quark mass near threshold. 
Our objective is that the definition should be gauge independent and well-defined within quan- 
tum field theory so that radiative and relativistic corrections can be calculated in a systematic 
way. Nevertheless, we stress that both concepts, as well as the other concepts mentioned before, 
are equivalently applicable to considerations near the quark production threshold. 

2.1 The PS mass 

We define 

mpg = mpoie — 5"^pg with Srrip^ = T.soit( 



(3) 



where Egoft is the soft part of the heavy quark self energy which is defined as the part where at 
least one of the heavy quark propagators is on-shell. A more precise definition is given below. 
We will show that the static PS mass and the new PS mass definition coincide at leading order 
in 1/m. At the same time we stress that the new definition accounts for recoil corrections of 
orders 1/m and But first we describe the physical picture and discuss the structure of 

the quark self energy shown in Fig. |I[ 

The starting point of our considerations is an on-shell quark with mass m and momentum p 
(i.e. J? = m^) which we consider to be at rest, p = (m, 0). This quark interacts with a number 
of gluons, the subdiagram 5* displayed in Fig. describes the interaction between the gluons. In 
general the quark lines between the interaction points represent virtual quark states. However, 
if the virtual quark comes very close to the mass shell and the total momentum of the cloud 
of virtual gluons becomes soft, this situation gives rise to long-distance nonperturbative QCD 
interactions. The described (virtual) contributions result in the soft part of the self energy, 
Ssoft- For a precise definition we start with a general self energy diagram as shown in Fig. 

- ^m) = / n 7^5g:|({/^}) (-^^^.t-^-t.^.o n ^--z Has^-r^j (4) 

m=l K'^'^J n=N rn '"' 
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where the last factor is a non-commutative product with decreasing index n. The hne momenta 
kn are hnear combinations of the gluon loop momenta Im, the particular representation is 
specified by the structure S. The symbol {Im} means the set of all these loop momenta, the 
same symbol is used for the Lorentz and color indices. In general we have M < N which means 
that line momenta can be correlated. The momenta of the virtual quark states are given by 
Pn = P + kn- Taking this as the starting point we define 

N M jAi i+l 

i=l"' m=l\'^'^> n=N Vn 

xz{^, + m)(-zn6{p^-m')){-zg,r'TaJ {[ {-igsT-T,^) ■ (5) 

^ ^ ntt-l Pn-m 

This equation is the definition of the soft part of the quark self energy. One can derive this 
expression from Eq. @) by using the identity 

1 _ ^ p / ^r^l (6) 



pZ _ jj^Z _|_ ypZ _ jy^. 

and the fact that the principal value integral does not give any infrared sensitive contribution. 
The delta function can be used to remove the integration over the zero component of fcj. In order 
to parameterize the softness of the gluon cloud we impose a cutoff on the spatial component, 
\ki\ < fif, and indicate this by a label fif written at the upper limit of the three-dimensional 
integral. This cutoff yU/ is also known as factorization scale. So we can rewrite Eq. (|^) as 



where 

^ M-l j4j i+l ■ 



m ' nt^LlPn-m 

The range of the index m is reduced by one which indicates that one of the loop momenta 
is extracted as line momentum of the i-th line. In the following we deal with the different 
realizations of this compact expression. As we will see explicitly, the function V{k,p) occurring 
as integrand can be seen as quark-antiquark potential where we have summed over the spin 
of the tensor product of a final state with an initial state. Because the static quark-antiquark 
potential is used in a similar way in Ref. we recover the result of Ref. [^] in the static 



limit. However, there is no kind of hierarchical order between both concepts because both solve 
the afore mentioned problems with non-perturbative uncertainties of the order O(Aqcd)- The 
role of non-perturbative effects of the order 1/m have to be studied elsewhere. 

2.2 The leading order perturbative contribution 

The leading order contribution to the self energy of the quark is given by 
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<x < ^ <x <^ <x 

(a) (b) (c) (d) 

Figure 2: Leading order contribution to the quark self energy (a) and to the quark- ant iquark 
potential (b). The cross indicates the point where we cut the quark line by imposing an on-shell 
condition to the virtual quark state. The gluon propagator can be decomposed in a Coulomb 
propagator (c) and a transverse propagator (d). 



where Feynman gauge is used for the gluon. The soft contribution thus reads 

= -^9lCp J (-2(^ + ^) + 4m) 6 {{p + k)' - m') . (10) 

The projector (1 + 7°)/2 which represents the on-shell quark at rest can be placed at the left 
and at the right of the Dirac structure and will lead to a further simplification. In terms of 
the zero component of the momentum k the Dirac delta function has two zeros ko = k^ and 
ko = k_ with 

k± := ibV + — m (11) 
where k = \k\. The delta function is therefore written as 

S ((p + ky - m') = J— ^ mo - k^) + 5{k, -k_)). (12) 



The procedure which is done here is shown in Fig. |^(a-b). The cross indicates that we cut the 
line at this point by imposing the on-shell condition to the corresponding (virtual) momentum. 
The diagram then proceeds to a quark-antiquark interaction diagram where we have kept the 
crosses to indicate the position of the cut line. This line carries the momentum p + k while the 
other two external lines carry the momentum p. Accordingly we obtain 



^soft I 



>/^/) = -2y (^^(^'^^)' V{Kp) = V^{k,p) + V.{k,p) (13) 



where 



ym^ + K'^{k^ — k"^) 2ni\/m? + K'^{\/m? + =]= m) 

For m <^ fif the restriction of the three-dimensional integral allows for an expansion in K/m. 
We obtain 

V.Ck,p) = -4na.cJ - + O (^)\ (15) 




and therefore 




Vikp) = -Ana^C^ <! - " ^ + O ( — I ^ (16) 



The first term is the Coulomb potential for a quark-antiquark interaction. The second term 
can be related to the Breit-Fermi potential of the quark-antiquark interaction by summing 
over the spin states of the tensor product of the final quark and the final antiquark state and 
using the same kinematic constraints. Moreover, one can identify V^(fc,p) with the scattering 
potential and V-{k,p) with the annihilation potential. 

One can integrate this radial symmetric potential over the space components and end up 
with an expansion in 1/m. Note that there are no contributions to odd powers of 1/m. But one 
can also integrate the exact expressions for V±{k,p) in Eq. (0) up to the factorization scale /i/ 
and obtain 

The expansion of this expression in small values of fif/m results in 

OsCf 



■J soft 



IT 




The first term reproduces the result given in Ref. to leading order in as while the second 
term is the recoil correction to the static limit in this order of perturbation theory. This second 
term is related to the Breit-Fermi potential but does not coincide with it. 



3 Two loop contributions 

To take a step beyond the leading order perturbation theory, we consider two-loop diagrams 
for the heavy quark self energy as shown in Fig. |^. We calculate them in Coulomb gauge, even 
though we stress that our final result is gauge invariant. The gluon propagator in Coulomb 
gauge is given by 

G^^{k) = '-^, ^?(^) = ^ (^^. - |t) ' ^'^■ = 1'2.3. (19) 
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Figure 4: Diagrams which cancel due to the classical Ward identity 
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(a) (b) (c) 

Figure 5: Diagrams which contribute to order 0(/ii/m^) 



The use of Coulomb gauge splits up the gluon propagators into a Coulomb term (Coulomb 
gluon) and a transverse term (transverse gluon) where the first one couples to the quark via 
the time components only. This splitting is shown in Fig. ^(b-d). 



3.1 The abelian diagrams 

We start our consideration with the abelian diagrams shown in Figs, ^(a) and (b). In cutting 
the quark line in all possible ways we obtain a lot of diagrams. However, we find that the 
final contribution of these diagrams to the soft part of the self energy is suppressed by fi'j/m'^. 
There can be found different arguments for this suppression. First, in applying the classical 
Ward identity, the QED diagrams shown in Fig. ^ cancel exactly at \k\ — 0, the remaining 
contribution is of order 0{k'^/m'^). Concerning this note that the Ward identity for the in- 
teraction vertex of a Coulomb gluon with the quark holds even in non-abelian theories P7|. 



A second argument is that the interaction between a transverse gluon and a non-relativistic 
quark as shown in Fig. |^(a) is suppressed by fif/m, leading to an overall fi'j/m'^ suppression. 
In addition, the box diagrams in Figs, ^(b) and (c) are either suppressed by a factor /i^/ 



m 



2 



or give an iteration of the leading order potential. To summarize, the diagrams in Figs. y(a) 
and (b) give contributions only of the order g^fi'j/m'^. 



3.2 The vacuum polarization of the gluon 

Following these arguments, it turns out that the only abelian diagrams which can give a non- 
suppressed contribution to the soft part of the quark self energy are the diagrams containing 
the vacuum polarization of the gluon as shown in Fig. |^(d-f). The simple calculation of these 
diagrams within the MS scheme, accounting only for light fermion loops, gluon loop (and ghost 
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(e) 

Figure 6: The non-abelian self energy diagram in Coulomb gauge; displayed are Coulomb gluons 
(wiggles) and transverse gluons (dashed lines) 



loop if Feynman gauge is used) results after renormalization in 



asifi) (sICa 20TfNf (HCa ATpNp 
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(3))) ■ 




W'|l + if'l''.-/3ol'M7llU- (20) 

This result has been anticipated because the expression in the curly brackets of the integrand 
reproduces the next-to-leading order correction to the QCD Coulomb potential. 

3.3 The non-abelian diagram 

In this subsection we calculate the non-abelian diagram shown in Fig. ^(c). In Coulomb gauge 
this diagram gives rise to seven two-loop diagrams which are shown in Fig. Direct calcula- 
tions show that only the diagram in Fig. ^(b) gives a contribution of order g'^fif/m while the 
other diagrams are of order g'^^'j/m'^ or vanish to this order in the coupling after applying the 
renormalization procedure (see e.g. Ref. |3^). The calculation of the diagram in Fig. |^(b) is 



simple and we show it in detail. The contribution of this diagram to the self energy is given by 

~ J {2nr kl kl {k, - k,f i^k, - k,f ) 

X (-^^?s7°T,)i±^. (21) 

2 f + f2 — rn p + fi — m 2 
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Figure 7: The soft part of the non-abehan diagram under consideration 



Here the Lorentz structure of the three- gluon vertex is reduced to ^1 + ^2- The Dirac structure 
of the integrand can be simphfied to 



1 + 7" 



{-(2m + k2oh%l - k2ll\2m + ho)} 



(22) 



The rest of the diagram is symmetric with respect to the interchange of the two hne momenta 
ki and /c2- That provides us with a further simphfication of the Dirac structure, the self energy 
reduces to 



g'^CpCA 



d'^ki d'^k2 



(4m + fcio + k2o) {klkl - (^1^2 



(27r)4 (27r)4 klkl{h - k2)\ki - k2)''{{p + h)^ - m^){{p + ^2)2 - 



where TaTbTcfabc = iCFCA/2 has been used. We now employ the replacements 
1 . . 1 



(23) 



{p + /ci)2 — m^ 



—in6{{p + ki) — m), 



{p + k2)^ — m? 



-m6{{p + k2)' - m') (24) 



i.e. we set cuts at the two intermediate quark lines separately to obtain the two parts ^^0^ 
and SgQ£^.2 of the soft contribution T^^^a '^^ self energy, as shown in Fig. 0. Taking the first 
cut, the delta function removes the integration over kiQ. At the same time our definition of 
the soft contribution imposes a restriction \ki\ < fif on the space components of the first line 
momentum. The delta function 



m 



1 



2Jm'^ + kl 



i5{kw - ki+) + 5{kw - ki_)) 



with 



ki± = ± Y + kj — m 



(25) 



(26) 
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causes two contributions which are known as scattering and annihilation amphtude (according 
to and ki-, resp.). The integration over /c2o is done by using the residue theorem. Actually 
there are only two denominator factors which can contribute to poles of the integrand, namely 
{ki — ^2)^ and ((p + ^2)^ ~ m^). In summing up the four contributions from the integration over 
k2o and kiQ we end up with a two-fold three-dimensional integral over the space components of 
the two line momenta where the first integral is restricted by \ki\ < fif as mentioned earlier. 
We now impose the restriction \ki\ < fif on the integrand to simplify it and obtain 

.66 _ gtCpCA rf d^k, f d^k2 klkl - ihh 



4m J {2-^)^] {27: fkl{kl?{h-k2Y' 

The integral over the space components of /c2 can be easily done by executing the angular 
integration followed by the radial integration. We obtain 



rf3fc2 klkl - {hhf 



(28) 



{2ny kl{k^f{ki-k2y I6\ki\ 
and therefore, finally 

y6b _ O^CfCa 2 (c.c^^ 



Symmetry considerations show that ^^^£^2 gives exactly the same contribution. As mentioned 
before, there are no other non-abelian contributions, therefore we obtain 

^Zi = ^^f^}. (30) 
This result has been anticipated, too, to be minus one half of the non-abelian correction to the 



QCD Coulomb potential, which is known in the literature (see for example Refs. |5 



.ATA 1 d^f^ i TT^af CiT'CA 1 a'^CpCA 2 



This calculation concludes the considerations of the two-loop diagrams shown in Fig. ^. 
3.4 Our final result 

Summarizing all contribution up to NNLO accuracy, we obtain 



TT I m Att \ mj \ An J 

(32) 

where m is the pole mass, /i is the renormalization scale. This result is considered to be of the 
order O(a^) because we imply that the ratio /i//m is typically of the order of or smaller. 
The scale ;U/ is the factorization scale, and 

Co = 1, C'q = 0, Cq = —-, 




C2 = a2-2{2a,Po + f3i)i\ni^yiyAp',[\n'l^^^ (33) 
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The constants ai, 02, Po, and jSi are given in Appendix A. The coefficients Ci and C2 have been 
derived in Ref. [|29l by using known corrections to the QCD potential. In this work we have 
derived the coefficients C'q, Cq, and C[. Note that our result can be represented in a condensed 
form as 

m^M-m = ~-J j^,{Vc{\k\) + Vn{\k\) + VM\k\)) (34) 

where the ffist term Vc is the static Coulomb potential, Vr is the relativistic correction (which 
is related to Breit-Fermi potential but does not coincide with it), and Vna is the non-abelian 
correction. 



4 Top quark production near threshold 
with NNLO accuracy 

In this section we consider the cross section of the process e^e^ —>■ ti in the near threshold 
region where the velocity v of the top quark is small. It is well-known that the conventional 
perturbative expansion does not work in the non-relativistic region because of the presence of 
the Coulomb singularities at small velocities f — 0. The terms proportional to (ag/v)"' appear 
due to the instantaneous Coulomb interaction between the top and the antitop quark. The 
standard technique for resumming the terms (a^/f)" consists in using the Schrodinger equa- 
tion for the Coulomb potential and to find the Green function 0, |^ . The Green function is then 
related to the cross section by the optical theorem. In order to determine NLO corrections to 
the cross section we need to know the short- distance correction to the vector current [H2l, the 



NLO correction to the Coulomb potential and the contribution of the non-factorizable cor- 



rections |jT5|, |T6|, |T^. It was proven that the non-factorizable corrections cancel in the inclusive 
cross section at NLO and beyond |T^, |16|, p!7| . 

At NNLO the situation is more complicated. One of the obstacles for a straightforward cal- 
culation are the UV divergences coming from relativistic corrections to the Coulomb potential. 
This problem can be solved by a proper factorization of the amplitudes and by employing effec- 
tive theories. We want to sketch the derivation of the inclusive cross section. The inclusive cross 
section can be obtained from the correlation function of two vector currents j^{x) = i[x)^^t{x), 



The ffist step is to expand the Lagrangian and the currents 



C2 



(35) 



(36) 



consistently in 1/m. The useful language for treating the one and two non-relativistic quark(s) 
is provided by the NRQCD ^ and the PNRQCD respectively. After the expansion, 
the relative cross section reads 



R = cr(e+e" ti)/apt = e|A^c— C(ro)Im 



1 



p' 



?>rn? 



G{ro,ro\E + iT) 



(37) 



ro^O 



where apt = 4:7ia'^/3s, eg is the electric charge of the top quark, Nc is the number of colors, 
y/s = 2m + E is the total center-of-mass energy of the quark-antiquark system, m is the top 
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quark pole mass and F is the top quark width. The unknown coefficient C(ro) can be fixed 
by using a direct QCD calculation of the vector vertex at NNLO in the so-called intermediate 
47|, H3] and by using the direct matching procedure suggested in Ref. [E^]. The result 



region 

of such a matching procedure for the coefficient C(ro) is given in the Appendix [IS, 20, 21, 22 



The function G{f, f'\E + iT) is the non-relativistic Green function. It satisfies the Schrodin- 
ger equation 

(H - E - iT)G{f,f'\E + iV) = 5{f-f'). (38) 
The non-relativistic Hamiltonian of the heavy quark-antiquark system reads 



H = Ho + W{r), Ho = ^ + Vcir), (39) 

m 

where Vc{r) is the static QCD potential of the heavy quark-antiquark system at NNLO order. 



Voir) 



+ 



r 



l + ^^^(2/?oln(/iV)+aO 



47r 



/?oM41n^(/iV 



TT 



2(/5i + 2/3oai)ln(/iV) + a2 



(40) 



yu' = ^e""^ 1 /i is the renormalization scale, and 7^; is Euler's constant. The color factors are 
given by Cp = 4/3, Ca = 3 and Tp = 1/2. The number of light quark fiavors is Np = 5. The 
coefficients Oi and 02 were calculated in and respectively, and are listed in the 

Appendix. 

The function W{r) in Eq. ( p9D is the QCD generalization of the QED Breit-Fermi Hamilto- 
nian I^G], |53|. We consider here the quark-antiquark production in the S'-wave mode. The 
Breit-Fermi Hamiltonian for the final state with L = 0, 5*^ = 2 reads 



W{r) 



P 



Am? 



+ 



\\-nCFCts{^i) 



6{r) 



Cpasifi) r 1 



2m^ 



r 



(41) 



It has been demonstrated in Refs. |20|, ^ that the Schrodinger equation can be reduced to 
the equation 

(Hi - Ei)Gi{r, r'l^i) = 5%r - r') 



with the energy Ei = E + E'^/Am, E = E + iT, and with the Hamiltonian 



H, = ^ + Vc{r) + ^Vo{r) 
m 2m 



3 ^ 



2m 



Voir) 



(42) 



(43) 



The final expression for the NNLO cross section is given by 



R 



^^i^o^E) = %\l-ACp^ 



Im 



TT 



bm I 



(44) 

with G'i(ro, rol-Ei) being the solution of Eq. (^2]) at r = r' = tq while C2{tq) is taken from 
Eqs. (|A^- |A5|) . For the numerical solution we use the program derived in Ref. by one of 
the authors. 
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m{m) 


mpg 


m^^o 


m^I'LO 

™NNLO 
AAipg 


'""PS 


"Sole 


™NLO 
''''pole 


™,NNLO 
'''pole 


00 

"^pole 


160.0 


166.51 


167.69 


167.97 
167.95 


168.05 
168.03 


167.44 


169.05 


169.56 


169.80 


165.0 


171.72 


172.93 


173.22 
173.20 


173.30 
173.28 


172.64 


174.28 


174.80 


175.05 


170.0 


176.92 


178.17 


178.47 
178.45 


178.55 
178.53 


177.84 


179.52 


180.05 


180.30 



Table 1: Top quark mass relations for the MS, PS, PS, and the pole mass at LO, NLO, NNLO, 
and "large /Sq" accuracy in GeV. We fix the MS mass to be m(m) = 160 GeV, 165 GeV, and 
170 GeV and find the pole mass at LO, NLO, and NNLO from the three-loop relation in 
Eq. (^Sf ). The PS and PS masses are derived from the pole mass by using Eq. (H) (without or 
with the 1/m contributions, respectively). We use the QCD coupling constant ^^(m^) = 0.119 
[p6| , /i = m(m), and fif = 20 GeV for the factorization scale in the PS and PS masses. 



4.1 Relations between the masses 



Our main result is Eq. (|32D, it enables us to relate the pole mass to the PS mass or the PS 
mass (the latter by leaving out the terms of higher order in 1/m). Note that we did not include 
electroweak corrections neither in this mass relation nor in the cross section. To relate the pole 



mass to the MS mass we use the three- loop relation which has been derived in Ref. |54, |55| 



""P"^^ - l + + (-1.0414iV^ + 13.4434) 



m(m) 3 V TT / V 7r 

+ f — y (0.6527A^| - 26.655iVF + 190.595) (45) 



TT 



where as = a;s(m) is taken at the MS mass. More precisely, we fix the MS mass to take the 
values m(m) = 160 GeV, 165 GeV, and 170 GeV and use Eq. ( ^3D to determine the pole mass 
at LO, NLO, and NNLO. This pole mass is then used as input parameter m in Eq. ( ^2]) to 
determine the PS and PS masses at LO, NLO, and NNLO. The obtained values for the PS and 
the PS mass differ only in NNLO. The results of these calculations are collected in Table |l|. 



together with the estimates for "large /3o" accuracy [^0|, ^ . Note that the same values for the 
MS mass have been used for Tables 2 and 3 in Ref. [26|. The obtained mass values can now be 
used for an analysis of the Schrodinger equation. 

Taking the PS mass instead of the PS mass, we observe an improvement of the convergence. 
The differences for the mass values in going from LO to NLO to NNLO to the "large /?o" 
estimate for e.g. m(m) = 165 GeV read 7.64 GeV, 1.64 GeV, 0.52 GeV, and 0.25 GeV for the 
pole mass, 6.72 GeV, 1.21 GeV, 0.29 GeV, anci 0.08 GeV for the PS mass and finally 6.72 GeV, 
1.21 GeV, 0.27 GeV, and 0.08 GeV for the PS mass^ 

In Fig. H we show the difference between the PS and the PS mass (in GeV) as a function 
of the factorization scale fij (solid line) at = 15 GeV. It is interesting to observe that the 
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non-abelian part of the difference between the PS and the PS mass (dotted hne in Fig. H) 
can be as large as 200 MeV. But the recoil correction cancel in part the non-abelian one and 
therefore the final difference is not more than 50 MeV. The dependence of mpg — mps on the 
renormalization scale at Hf = 30 GeV is given in Fig. ^ by the dashed curve. 



4.2 The concept with the pole mass 

The top quark cross section at LO, NLO, and NNLO is shown in Fig. ^ as a function of the 
center-of-mass energy. For the top quark pole mass we use rrit = 175.05 GeV, for the top quark 



width Tt = 1.43 GeV, and for the QCD coupling constant as{mz) = 0.119 [§6l. Different 
values fi = 15 GeV, 30 GeV, and 60 GeV for the renormalization scale are selected because 
they roughly correspond to the typical spatial momenta for the top quark. For solving the 
Schrodinger equation we use the program written by one of the authors ||2l|. Note that we do 
not take into account an initial photon radiation which would change the shape of the cross 
section. This can be easily included in the Monte Carlo simulation. 

The NNLO curve modifies the line shape by the amount of 20 — 30% which is as large as the 
NLO correction. It also shifts the position of the 15* peak by approximately 600 MeV. These 
large shifts of the peak position were expected. As we discussed above (and is well-known in 



the literature [^, |28|), the pole mass definition suffers from the renormalon ambiguity. The 
top quark pole mass cannot be defined better than O(Aqcd)- Large NNLO corrections and a 
large shift of the IS" resonance can spoil the top quark mass measurement at the NLC. 



4.3 The concept with the PS mass 

In this subsection we discuss the calculation concept using the PS mass. We redefine the pole 
mass through the PS mass using the relation given in Eq. (|32D without the l/m contributions 
and then use the PS mass as an input parameter for our numerical analysis at LO, NLO, and 
NNLO. In Fig. |10| we show the top quark cross section expressed in terms of mps(/i/) at LO, 
NLO, and NNLO (like in Fig. |^) as a function of the center-of-mass energy. We take mps(/i/ = 
20 GeV) = 173.30 GeV which corresponds to the pole mass m = 175.05 GeV according to 
Table |1|. In looking at Fig. |10] we observe an improvement in the stability of the position of the 
first peak in comparison to the previous analysis as we go from LO to NLO to NNLO. Actually, 
for the three values /i = 15 GeV, 30 GeV, and 60 GeV we obtain the maxima of the IS* peak for 
NNLO at Sinax = 347.32 GeV, 347.41 GeV, and 347.48 GeV while the maximal values are given 
by -Rmax = 1.379, 1.184, and 1.088, respectively. This demonstrated that a large variation in 
the renormalization scale fi gives rise only to a shift of about 160 MeV for the IS* peak position 
at NNLO while the variation for -Rmax is still large. 



4.4 The concept with the PS mass 

In this subsection we discuss the calculation concept where we use the PS mass. We redefine 
the pole mass through the PS mass by using the relation given in Eq. (|32|) and then use the PS 
mass as an input parameter for the numerical analysis at LO, NLO, and NNLO. In Fig. ^ we 
show the top quark cross section expressed in terms of mpg{nf) at LO, NLO, and NNLO (like 
in Fig. 1^) as a function of the center-of-mass energy. We take ?7T.pg(/i/ = 20 GeV) = 173.28 GeV. 
Again we observe a very good stability of the position of the first peak as we go from LO to 
NLO to NNLO, similar to the one observed for the PS mass case. Studying the size of the 
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NNLO corrections to the peak positions we conclude that the current theoretical uncertainty 
of the determination of the PS mass from the 15* peak position is about 100 — 200 MeV. 



5 Conclusion and discussions 

We introduced the so-called PS mass which is an infrared-insensitive quark mass, naturally de- 
fined by subtracting the soft part of the quark self energy from the pole mass. We demonstrated 
a deep relation of this definition with the static quark- ant iquark potential. At leading order 
in 1/m the PS mass coincides with the PS mass, defined in a completely different manner. In 
contrast to the pole mass, the PS mass is not sensitive to the non-perturbative QCD effects. 
Going beyond the static limit, the small normalization point introduces recoil corrections to 
the relation of the pole mass with the PS mass which have been calculated in this paper. In 
addition we demonstrated that, if we use the PS or the PS mass in the calculations, the per- 
turbative predictions for the cross section become much more stable at higher orders of QCD 
(shifts are below 0.1 GeV). This understanding removes one of the obstacles for an accurate 
top mass measurement and one can expect that the top quark mass will be extracted from a 
threshold scan at NLC with an accuracy of about 100 — 200 MeV. 
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Appendix 



The coefficients ai and 02 were calculated in Refs. and |]^, ^ , respectively, and are given 
by 

31 20 

/4343 , 2 22 \ ^2 56 

02 = + 47r2 + — Cs ]Cl- [ + — Cs CaTfNf 

V 162 4 3 ^7 ^ V 81 3 ^7 ^ ^ ^ 

20_ _ \2 /55 



+ [-^TpNpj - (^y - leCsj CfTfNf. (Al) 

The first two coefficients in the expansion of the QCD /5-function are 

11 4 34 20 

/3o = ^Ca - -TfNf, I3i = —C\ - —CaTfNf - ACfTfNf. (A2) 

The short distance coefficient reads 

dr) = l-ACF^ + C2ir)CF(^]\ 

71 y TT J 

C2{r) = Ai\og{r/a) + A2\og{m/fih) + As (A3) 
where the hard renormalization scale is taken to be the pole mass, Hh = m, a = e^~""^ /2m, and 

Al = TT^ [Ca + ^) , A2 = 2/3o, 

A3 = C^Cf + C^^Ca + C^TfNf + C^Tf (A4) 
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where Np = 5 is the number of the hght flavors, and 

39 ^ 2/4, ^ 35 

= Cs + TT - n 2 

^ 4 V3 18 



151 13. 0/179 8 
36 2 ■ V 72 3 



= -^-^C3 + vrM^-^ln2 



= 11, = 11 - 1.^ (A5) 

2 9 ' 2 g g V ; 
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Figure 8: The difference between the PS and the PS mass (in GeV) as a function of the 
factorization scale fif (sohd hne) at /i = 15 GeV. The dotted hne shows only the non-abelian 
part of the difference. The dependence of m-p^ — mps as a function of the normalization scale 
/i at /ij = 30 GeV is shown as dashed line. 
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Figure 9: The concept using the pole mass: shown is the relative cross section R{e~^e~ — > tt) as 
a function of the center-of-mass energy in GeV for the LO (dashed-dotted lines), NLO (dashed 
lines), and NNLO (solid lines) approximation. We take the value mt = 175.05 GeV for the pole 
mass of the top quark, Tt = 1.43 GeV for the top quark width, as{mz) = 0.119 and different 
values /i = 15 GeV, 30 GeV, and 60 GeV for the renormalization scale. 
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Figure 10: The concept using the PS mass: shown is the relative cross section R{e~^e~ tt) as 
a function of the center-of-mass energy in GeV for the LO (dashed-dotted hues), NLO (dashed 
hues), and NNLO (sohd hues) approximation. We take the value mps = 173.30 GeV for the PS 
mass of the top quark, Tt = 1.43 GeV for the top quark width, as{mz) = 0.119 and different 
values fi = 15 GeV, 30 GeV, and 60 GeV for the renormalization scale. 
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Figure 11: The concept using the PS mass: shown is the relative cross section R{e^e^ tt) as 
a function of the center-of-mass energy in GeV for the LO (dashed-dotted lines), NLO (dashed 
lines), and NNLO (solid lines) approximation. We take the value frip^ = 173.28 GeV for the PS 
mass of the top quark, P^ = 1.43 GeV for the top quark width, as{mz) = 0.119 and different 
values /i = 15 GeV, 30 GeV, and 60 GeV for the renormalization scale. 
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